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Abstract 

This work is focused on a nonlinear equation describing the oscil- 
lations of an extensible viscoelastic beam with fixed ends, subject to 
distributed elastic external force. For a general axial load j3, the exis- 
tence of a finite/infinite set of stationary solutions and buckling occur- 
rence are scrutinized. The exponential stability of the straight position 
is discussed. Finally, the related dynamical system in the history space 
framework is shown to possess a regular global attractor. 
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1 Introduction 



In this paper we analyze the asymptotic behavior of the following nonlinear 
dissipative evolution problem 

/•oo 

d t tu+d xxxx u+ / fi(s)d xxxx [u(t)-u(t-s)}ds-((3+\\d x u\\ 2 L 2^ _ 1) )d xx u = -ku+f, 
Jo 

(1) 

in the unknown variable u — u(x,t) : [0, 1] X K. — > K, which describes the 
vibrations of an extensible viscoelastic beam of unitary natural length. The 
real function / = f(x) is the lateral (static) load distribution, the term — ku 
represents the lateral action effected by the elastic foundation and the real pa- 
rameter j3 denotes the axial force acting in the reference configuration (positive 
in stretching, negative in compression) . The memory kernel [L is a nonnegative 
absolutely continuous function on ]R + = (0,oo) (hence, differentiable almost 
everywhere) such that 

p!(s) + 6fi(s) < 0, / (i(s)ds = k, (2) 
Jo 

for some S > and n > 0. In particular, /x has an exponential decay of rate S 
at infinity. 

Precisely, the model is obtained by combining the pioneering ideas of Woinowsky- 
Krieger [53] with the theory of Drozdov and Kolmanovskii [TJ], i.e. taking into 
account the geometric nonlinearity due to the deflection, which produces an 
elongation of the bar, and the energy loss due to the internal dissipation of the 
material, which translates into a linear viscoelastic response in bending. 

In our analysis, we assume that both ends of the beam are hinged; namely, 
for every t E M. 

u(0, t) = u(l, t) = d xx u{0, t) = d xx u{l,t) = 0. (3) 

Moreover, because of integro-differential nature of (p}, the past history of u 
(which need not fulfill the equation for negative times) is assumed to be known. 
Hence, the initial condition reads 

u(x,t) — v,i,{x,t), (x,t) € [0, 1] x (-oo,0], (4) 

where u+ : [0, 1] x (— oo, 0] — > K is a given function. 

In order to apply the theory of strongly continuous semigroups, we recast 
the original problem as a differential system in the history space framework. 
To this end, following Dafermos we introduce the relative displacement 
history 

rf(x, s) — u(x, t) — u(x, t — s), 



so that equation (JT|) turns into 



d tt u + d xxxx u + / ^(s)d xxxx n(s)ds - (f3 + \\d x u\\ L2 / u)d xx u 
Jo 

^d t r) = -d s n + d t u. 

Accordingly, the initial condition (0} becomes 

u(x, 0) = u (x), i€ [0,1], 
d t u(x, 0) = u%(x), i£ [0,1], 
k rf (x, s) = r) (x, s), (x, s) e [0, 1] x R+ , 

where we set 



-ku + /, 
(5) 

(6) 



uq(x) = u±(x, 0), = d t Ui,{x, 0), s) = u*(a;, 0) - -s). 

As far as the boundary conditions are concerned, ([3]), for every t > 0, translates 
into 

f «(0,t) = U(l,t) = flSrxUfO,*) = fllrxUfl,*) = 0, 

fj*(0, s) - fftl, s) = a M r/*(0, s) = 9 ra r/*(l, s) = 0, (7) 

rf(x, 0) = lim n (x, s) = 0. 

s— >o 

It is worth noting that the static counterpart of problem (TTJ) reduces to 



W 



P+\\d x u\\ 



u 



L 2 (0,1) 

u(0) = u(l) = u"(0) - u"(l) = 



ku = /, 



(8) 



When k = the investigation of the solutions to (O and their stability, in 
dependence on /3, represents a classical nonlinear buckling problem (see, for 
instance, [U [T3J [50]). Its numerical solutions are available in [7J, while their 
connection with industrial applications is discussed in 15!. Recently, a careful 
analysis of the corresponding buckled stationary states was performed in [19] 
for all values of /?. In [9] this analysis was improved to include a source / with 
a general shape. 

For every k > and vanishing sources, exact solutions to ((8]) can be found 
in [B], and at a first sight, this case looks like a slight modification of previ- 
ously scrutinized models where k vanishes. This is partially true. Indeed, the 
restoring elastic force, acting on each point of the beam, opposes the buckling 
phenomenon and increases the critical Euler buckling value j3 c , which is no 
longer equal to y/Xi, the root square of the first eingenvalue of the d xxxx oper- 
ator, but turns out to be a piecewise linear function of k. When the lateral load 
/ vanishes, the null solution is unique provided that /3 > —fi c (k) (see Theorem 
[2]), and buckles when (3 exceeds this critical value, as well as in the case k = 0. 



On the contrary, for some special positive values of k, called resonant values, 
infinitely many solutions may occur. 

Moreover, in the case k = 0, if / vanishes, the exponential decay of the en- 
ergy is provided when j3 > — Ai , so that the unique null solution is exponentially 
stable. On the contrary, as the axial load j3 < — Ai the straight position loses 
stability and the beam buckles. So, a finite number of buckled solutions occurs 
and the global attractor coincides with the unstable trajectories connecting 
them. 

By paralleling the results for k = 0, the null solution is expected to be 
exponentially stable, when it is unique. Quite surprisingly, it is not so. For large 
values of k, the energy decays with a sub-exponential rate when — (3 > (3 > — j3 c 
(see Theorem[8]). In particular, for any fixed k > Ai, the positive limiting value 
/3(k) is smaller than the critical value P c {k), and the first overlaps the latter 
only if < fc < Ai. 

The motion equation ((T|) with fi = k = / = turns out to be conservative 
and has been considered for hinged ends in ^ I12j , with particular reference to 
well-posedness results. Adding an external viscous damping term Sdtu (5 > 0) 
to this conservative model, stability properties of the unbuckled (trivial) states 
have been established in [3] U3] and, more formally, in [5D]. In this case, the 
global dynamics of solutions for a general j3 has been tackled first in [T5] . 
where some regularity of the attractor is obtained provided that S is large 
enough. When the stiffness of the surrounding medium is neglected (k = 0), the 
existence of a regular attractor was proved for extensible Kirchhoff beam [10] , 
extensible viscoelastic [17] and thermo-elastic [16] beams. A similar result for 
an extensible elastic beam resting on a viscoelastic foundation was obtained in 
[6] . This strategy can be generalized to the investigation of nonlinear dissipative 
models which describe the vibrations of extensible thermoelastic beams [3] , and 
to the analysis of the long term damped dynamics of extensible elastic bridges 
suspended by flexible and elastic cables In the last case the term — ku 
is replaced by —ku + and it represents a restoring force due to the cables. 
Moreover, our approach may be adapted to the study of simply supported 
bridges subjected to moving vertical load [32]. 

The final result of this work concern the existence of a regular global at- 
tractor for all values of the real parameter /3. The main difficulty comes from 
the very weak dissipation exhibited by the model, entirely contributed by the 
memory term. So, the existence of the global attractor is stated through the 
existence of a Lyapunov functional and the asymptotic smoothing property of 
the semigroup generated by the abstract problem via a suitable decomposition 
first proposed in [FT] . 

2 The Abstract Setting 

In this section we will consider an abstract version of problem ([S])-®. To 



this aim, let (Ho, (•, •), || -||) be a real Hilbert space, and let A : T>(A) <s H — » H a 
be a strictly positive selfadjoint operator. For £ £ R, we introduce the scale of 
Hilbert spaces 

fli = 23(^/ 4 ), (u,v) t ={A t /\A t /*v), \\u\U = \\A^u\\. 
In particular, Hg + i g iZ^ and the generalized Poincare inequalities hold 

>AilHI?<IHI?+ii Vu£H e+ll (9) 

where Ai > is the first eigenvalue of A. 

Given n satisfying we consider the L 2 -weighted spaces 

/>00 

^ = L*(]R+,^ +2 ), (r),Ot,„ = / Ai(*)M«)»£(«))<+2<k, = (r,,r])^ 

Jo 

along with the infinitesimal generator of the right-translation semigroup on 
A^Oi that is, the linear operator 

Tr] = -Drj, V(T) = {q £ M : Drj £ Mo, ??(0) = 0}, 

where D stands for the distributional derivative, and ?y(0) = lim^o r)(s) in H 2 . 
Besides, we denote by M\ the weighted Sobolev spaces 

M\ = Hl(R + ,H e+2 ) = { v £Mr-D V £ Me}, \\r,f M i = hll^ + \\D V \\l^. 

Moreover, the functional 

J(v) = - n'WWvWWld* 



v\\in ^ ( 10 ) 



is finite provided that rj £ T>(T). From the assumption (U) on fi 

1 

5 

Finally, we define the product Hilbert spaces 

Hi = H e+2 x Hi x Mi. 

For f3 £ M and / £ Hq, we investigate the evolution system on Hq in the 
unknowns u(t) : [0, oo) -> H 2 , d t u(t) : [0, oo) -» ifo and r/* : [0, oo) -» A^o 

' poo 

d tt u + Au+ I n^A^ds + (f3 + \\u\\f)A 1/2 u = -ku + f, 
d t r] = Trj + d t u, 
with initial conditions 

(li(0),U t (0),?7 ) = (u ,Ui,1] ) = z £ Hq. 



Remark 1. Problem (0-([7|) is just a particular case of the abstract system 
(fTTj) , obtained by setting H a = L 2 (0, 1) and 

A = d xxxx , V{d xxxx ) = {w € H 4 (0, 1) : w(0) = w(l) = w"(0) = w"(l) = 0} 

This operator is strictly positive selfadjoint with compact inverse, and its dis- 
crete spectrum is given by A n = n 4 7r 4 , neN. Thus, Ai = 7r 4 is the smallest 
eigenvalue. Besides, the peculiar relation [d xxxx ) x ^ 2 = —d xx holds true, with 

v(-d xx ) = h 2 (o, i)nfl£(o,i). 

Besides, system (jlip generates a strongly continuous semigroup (or dynam- 
ical system) S(t) on Ho which continuously depends on the initial data: for 
any initial data z € Jig, S{t)z is the unique weak solution to with related 
(twice the) energy given by 

£{t) - ||S(t)*||!i = + \\dtu(t)\\ 2 + \\r,%^. 

We omit the proof of these facts, which can be demonstrated either by means 
of a Galerkin procedure or with a standard fixed point method. In both cases, 
it is crucial to have uniform energy estimates on any finite time-interval. 



3 Steady states 

In the concrete problem ([5])-([7]), taking for simplicity / = 0, the stationary 
solutions (m, 0,0) solve the boundary value problem 



(12) 



U""- (/3+iKiii W) y'+te=o, 

\u(0) = u(l) = u"(0) = u"{l) = 0. 

Our aim is to analyze the multiplicity of such solutions. 
For every k > 0, let 

k 

Mn(fc) = —5—5- + n 2 tt 2 , P c (k) = min^„(fc). 

n 7T neN 

Assuming that nk € N be such that fi n = min/i„(fc), then it satisfies 

neN 

(n k -l) 2 n 2 < ^ <n 2 (n k + l) 2 . 

As a consequence, j3 c {k) is a piecewise-linear function of k. 
We consider the resonant set 

n = {i 2 j 2 Tr 4 :i,jeN,i <j}. 

When k G 1Z there exists at least a value fJ-j{k) which is not simple and \ii = 
i 7^ j, provided that k = i 2 j 2 7r 4 (resonant values). In the sequel, let /im(k) be 



the smallest value of {/x n } ra eN which is not simple. Of course, the /z„(fc) are all 
simple and increasingly ordered with n whenever k < 47r 4 . Given k > 0, let n* 
be the integer- valued function given by 

n^) = \Np\, Np = {n eN:/3 + M n(fc) < 0}, 

where |7V| stands for the cardinality of the set N . 

In the homogeneous case, we are able to establish the exact number of 
stationary solutions and their explicit form. In particular, we show that there is 
always at least one solution, and at most a finite number of solutions, whenever 
the values of fJ, n (k) not exceeding —fj are simple. 

Theorem 2. (see [%) If (i > —(3 c (k), then for every k > system (fT2)) has 
only the null solution, depicting the straight equilibrium position. Otherwise: 



if k £ K and /3 < — /i m (fc), the smallest non simple eigenvalue, there are 
infinitely many solutions; 

if either k £ TZ and -u. m (k) < (3 < —f3 c {k), or k ^ TZ and f3 < —p c (k), 
then besides the null solution there are also 2n*(/3) buckled solutions, 
namely 

u n( x ) — sin(n7ra) , n = 1, 2, . . . , n* (13) 



with 



1 



At = ± — y/-2 [(3 + /i n (k)]. (14) 



n 7r 



Remark 3. Assuming k — we recover the results of JIJ. 

When k ^ TZ the set of all stationary states is finite. Depending on the value of 
k and (3, the solutions branch in the pairs from the unbuckled state A^ — at 
the critical value /3 = — /3 c (fc), i.e., the beam can buckle in either the positive or 
negative directions of the transverse displacement. These branches exist for all 
/3 < —fi c (k) and A^ are monotone increasing functions of \/3\. For each n, (fl4l 
admits real (buckled) solutions A^ if and only if (3 < — /i„. When k < 47r 4 , for 
any /? in the interval 

<n + l) 2 TT 2 < p < ^ - n 2 TT 2 



(n + 1) 2 ^ 2 



the set Sq of the stationary solutions contains exactly 2n^ + l stationary points: 
the null solution and the solutions represented by ([13]) . 

When k £ TZ the set Sq contains an infinite numbers of solutions and all 
the resonant values are obtained by solving 



m, n £ N , n > m . 



The smallest value k is then equal to 4ir 4 and occurs when n = 2, m = 1. In 
the sequel we present a sketch of different bifurcation pictures occurring when 
k = n 4 (see Fig.l), and k = 97r 4 (see Fig.2). 




Fig.l. The bifurcation picture when k = it 4 




Fig.2. The bifurcation picture when k — 9ty 4 



4 The Lyapunov functional 

It is well known that the absorbing set gives a first rough estimate of the 
dissipativity of the system. In addition, it is the preliminary step to scrutinize 



its asymptotic dynamics and hence to prove the existence of a global attractor. 
Unfortunately, when the dissipation is very weak, a direct proof via explicit 
energy estimates might be very hard to find. For a quite general class of the 
so-called gradient systems it is possible to use an alternative approach appealing 
to the existence of a Lyapunov functional. This technique has been adopted in 

Definition 4. The Lyapunov functional is a function C E C(Ho,M) satis- 
fying the following conditions: 

(i) C{z) — > +00 if and only if ||z||« — > +00; 

(ii) C(S(t)z) is nonincreasing for any z € Ho; 

(iii) C{S{t)z) = C(z) for all t > implies that zeS. 
Proposition 5. The function 

C(t) = £{t) + IK*)lli) 2 + k \\u(t)\\ 2 2 (/, u(t)) 

is a Lyapunov functional for S(t). 

Proof. The continuity of £ and assertion (i) above are clear. Using (fTT|) . 
we obtain quite directly the inequality 

jC{S{t)z) < -dUWl^ (15) 

which proves the monotonicity of C along the trajectories departing from z. 
Finally, if C(S(t)z) is constant in time, we have that rf = for all t, which 
implies that u(t) is constant. Hence, z = S(t)z = (u ,0,0) for all t, that is, 
z eS. □ 

The existence of a Lyapunov functional ensures that 

Lemma 6. For all t > and initial data z 6 Wo, with ||z||-h < R, there 
exists a positive constant C (depending on \\f\\ and R) such that 

£(t) < C. 

Proof. Inequality (fT5)l ensures that 

C(t)<C(0)<Co(R,\\f\\) ■ 
Moreover, taking into account that 

IK*)f<f IK*)ll2<f m = Cie{t) 



□ 



we obtain the estimate 

L{t) > E{t) -2 {f, «(*)) > E{t) \ H/ll 2 - e \\u(t)\\ 2 > (1 - e d)5(t) - J ||/|| 
Finally, fixing e < we have 

f (t) < + J H/ll 2 ) ^ tttcT ( Co(i? ' + \ ll/l12 ) = c ' 

Moreover, defined the following functional 

$(t) = £(t) + £ (ftu , U) , 

where E(t) = £{t) + §(/? + IK*)Hi) 2 + k \\u(t)\\ 2 , we can prove 

Lemma 7. For any given z G Ho and for any t > and (3 e M, when e is 
small enough there exist three positive constants, mo, mi and m,2, independent 
of t such that 

m £(t) < $(t) < mi£(f) + m 2 . 

Proof. In order to prove the lower inequality we must observe that, by 
Young inequality 

\{dtu,u)\>-\ M 2 -^NI 2 ; 

hence, we obtain 



it) > \\u( t )\\i + (i - 1) im*)ii 2 + \ (p + \\u(t)\\iy 

(*-§) \Ht)f + \\v%,. 



+ 



If we choose e small enough to satisfy e < 2 and £ < 2k, then we have 

$(*) > m £(i) > m S(t), (16) 



where mo = min{l — |, 1 — ^}. 

The upper inequality can be obtained using the definition of $ and applying 
the estimate 

First, we can write 



\{dtu,u)\<-\\d t u\\ 2 + — \\u\\l ■ (17) 



*(*)< 



+ TT ( fe + |)]H*)ii2 + (| + 1 )ii«)ir 

+ ^(/3+IKi)||?) 2 + ||r/*||^. 



Then, by © and Lemma [5] we infer 



( ( 8 + lK)<| ( 8| + - 7 =C = C, (18) 
so that we finally obtain 

£(*) + \c 2 = mi£(t)+m 2 . 

□ 

As a byproduct, we deduce the existence of a bounded absorbing set Bo, chosen 
to be the ball of Ho centered at zero of radius i?o = 1+sup { ||^||ko : £( z ) — 
where K = l+sup 2oG5 C(zq). Note that Rq can be explicitly calculated in terms 
of the structural quantities of our system. 



5 Exponential stability 



Recalling Theorem [2J the set So of stationary solutions reduces to a singleton 
when 



P > -fic(k) = - min^„(fc), ft n (k) = y/X 



fc 



A„ = nV. (19) 



It is worth noting that /3 c (fc) is a piecewise-linear function of fc, in that fi c {k) 
fii(k) when < fc < vAiV^2 = 47r 4 , and in general 



/3 c (fc) = /i„(fc) when ^Xn-iVK < k < \/Xn^/xT l 



i+i- 



Unlike the case fc = 0, the energy £(t) does not decay exponentially in the 
whole domain of the (/3, fc) plain where ()19|) is satisfied, but in a region which 
is strictly included in it. 
More precisely, let 



Pc{k), < fc < Ai, 

2y/k, fc > Ai, 

the following result holds 

Theorem 8. When f = 0, the solutions to ([1]) decay exponentially, i.e. 

£{t) < c o £(0)e- ct 
with cq and c suitable positive constants, if and only if j3 > —/3(k). 



Using the same strategy bolstered in [SJ [B], the proof of this theorem is a 
direct consequence of the following lemma: 

Lemma 9. (see [5]) Let (3 € R, k > and 

Lu = Au + (3A^u + ku . 

There exists a rea7 function v — v(fi, k) such that 

(Lu, u) > v\\u\\ 2 , 

where v((3, k) > if and only if f3 > —f3(k). 

Remark 10. We stress that Theorem holds even if k = 0. In this 
case however we have (3(0) = — /3 C (0) = — yAi. Then the null solution is 
exponentially stable, if unique. 



6 The Global Attractor 

Now we state the existence of a global attractor for S(t), for any fj E R and 
k > 0. We recall that the global attractor A is the unique compact subset of 
Hq which is at the same time, fully invariant, i.e. S(t)A = A, for every t > 
and attracting, i.e. 

Urn S(S(t)B,A) ->0, 

t— >oo 

for every bounded set B C Ho, where 6 stands for the Hausdorff semidistance 
in Ho (see [J, [H], HQ)- 

We shall prove the following 

Theorem 11. The semigroup S(t) on Ho possesses a connected global at- 
tractor A bounded in H2, whose third component is included in T>(T), bounded 
in A4q and pointwise bounded in H4. Moreover, A coincides with the unstable 
manifold of the set S of the stationary points of S(t), namely, 

{_ z is a complete (bounded) trajectory of S(t) : 1 

gnM-t)-S\\ Ho = ]■ 

The set S of all stationary solutions consists of the vectors of the form (it, 0, 0), 
where u is a (weak) solution to the equation 

Au + (/? + \\u\\ 2 ^)A 1/2 u + ku = f. 

It is then apparent that S is bounded in Ho- If S is finite, then 

A = {2(0) : lim \\S(-t) - Zl \\n a = lim \\z(t) - z 2 \\ Ha - o}, (20) 



for some Z\, z% G S. 



Remark 12. When f = 0, k > and (3 > -f3 c (k), then A = S = 
{(0,0,0)}. If (3 < —j3 c (k), then S = So may be finite or infinite, according to 
Theorem [5) In the former case, (|2T)]) applies. 

The existence of a Lyapunov functional, along with the fact that S is a 
bounded set, allow us prove the existence of the attractor exploiting a general 
result from [8], tailored for our particular case. 

Lemma 13. (see [8]) Assume that, for every R > 0, there exist a positive 
function ipR vanishing at infinity and a compact set JCr C Hq such that the 
semigroup S(t) can be split into the sum L(t) + K(t), where the one-parameter 
operators Lit) and K(t) fulhll 

\\L{t)z\\ Uo <^R{t) and K{t)z G Kr, 

whenever \\z\\u < R and t > 0. Then, S(t) possesses a connected global 
attractor A, which consists of the unstable manifold of the set S. 

The proof of Theorem QT] will be carried out be showing a suitable asymp- 
totic compactness property of the semigroup, obtained exploiting a particular 
decomposition of S(t) devised in [T7] . 
By the interpolation inequality 

N|?<HH| 2 , 

it is clear that, provided that 7 > is large enough, 

~IM|f < \\u\\l+P\Hl+lM\ 2 <m\\u\\l (21) 

for some m — m(/3, 7) > 1. Now, choosing 7 = a + k, where k > is a fixed 
value, we assume a large enough so that (|2"T|) holds true. Let R > be fixed 
and || z || "H < R- Paralleling the procedure given in [17] . we decompose the 
solution S(t)z into the sum 

S(t)z = L(t)z + K(t)z, 

where 

L(t)z = (v(t),d t v(t),e) and K(t)z = (w(t),d t w(t),C) 
solve the systems 

{d tt v + Av + J Li(s)M(s)ds + (0 + \\u\\ \)A 1/2 v + av + kv = 0, 
d t i = T£ + d°v, ( 22 ) 
(v(0),d t v(Q),e)=z, 



and 



{d tt w + Aw + J n(s)A((s)ds + (/? + ||u||?)A 1/2 u; - av + kw = f, 
d t C = TC + dtw, ( 23 ) 
KO),^(0),C°) = 0. 

Then, Theorem [11] is proved as a consequence of the following lemmas. 

Lemma 14. There is lu = to(R) > such that 

\\L(t)z\\ no <Ce-^. 

It shows the exponential decay of L(t)z by means of a dissipation integral (see 
Lemma 5.2 of [17]). 

Lemma 15. The estimate 

\\K(t)z\\ Ha <C 

holds for every t > 0. 

It shows the asymptotic smoothing property of K(t) in a more regular space, 
for initial data bounded by R (see Lemma 6.3 of [17]). 

The proof of these two lemmas can be done following the guidelines of [T7] , 
remembering that in definition of the coefficient a in the functional $o (see 
[T7] pag. 726) also the elastic coefficient k must be taken into account. 
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